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ABSTRACT 



We report on our analysis of XMM-Newton observations of the Seyfert 2 galaxy ESO 138-Gl(z = 0.0091). These data reveal a 
complex spectrum in both its soft and hard portions. The 0.5-2 keV band is characterized by a strong "soft-excess" component with 
several emission lines, as commonly observed in other narrow-line AGN. Above 3 keV, a power-law fit yields a very flat slope (F ~ 
0.35), along with the presence of a prominent line-like emission feature around ~6.4 keV. This indicates heavy obscuration along the 
line of sight to the nucleus. We find an excellent fit to the 3-10 keV continuum with a pure reflection model, which provides strong 
evidence of a Compton-thick screen, preventing direct detection of the intrinsic nuclear X-ray emission. Although a model consisting 
of a power law transmitted through an absorber with Nh ~ 2.5 x 10"' cm"- also provides a reasonable fit to the hard X-ray data, the 
equivalent width (EW) value of ~ 800 eV measured for the Fe Ko- emission line is inconsistent with a primary continuum obscured 
by a Compton-thin column density. Furthermore, the ratio of 2-10 keV to de-reddened [OIII] fluxes for ESO 138-Gl agrees with the 
typical values reported for well-studied Compton-thick Seyfert galaxies. Finally, we also note that the upper limits to the 15-150 keV 
flux provided by Swi/z/BAT and INTEGRAL/lBlS seem to rule out the presence of a transmitted component of the nuclear continuum 
even in this very hard X-ray band, hence imply that the column density of the absorber could be as high as 10^^ cm This makes 
ESO 138-Gl a very interesting, heavy Compton-thick AGN candidate for the next X-ray missions with spectroscopic and imaging 
capabilities above 10 keV. 

Key words. Galaxies: individual: ESO 138-Gl - Galaxies: active - Galaxies: nuclei - X-ray: galaxies 
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Since X-ray radiation is emitted from within a few gravitational 
radii of the SMBH, it can provide insights into AGN activity 
with unrivalled detail. Furthermore, typical X-ray observations 
in the 2-10 keV range easily penetrate regions (i.e. dust lanes, 
the torus, star-forming knots) with high extinction in optical (Ay 
> 50 mag) located along our line-of-sight to the nucleus. The 
census of AGN in the Universe is, however, still incomplete as 
illustrated by the mismatch between the total mass function of 
local SMBHs and the SMBH mass function inferred from the 
X-ray-selected AGN luminosity functions (Marconi et al. 2004), 
and because only ~ 60% of the cosmic X-ray background in 
the 5-10 keV band has been so far resolved into point sources 
(Worsley et al. 2005). 

It is well-known that one of the major difficulties to over- 
come is the detection of the most elusive AGN, i.e. those suf- 
fering heavy obscuration (i.e. with column density of Nh i cry^ 
^1.6 X lO^"* cm-^), the so-called Compton-thick (CT) AGN. As 
the primary continuum is blocked, their X-ray emission below 
10 keV is mainly due to the Compton reflection and fluores- 
cence from optically-thick circumnuclear material (Matt 1997). 
Accordingly, the observed 2-10 keV fluxes for CT AGN there- 
fore represent only 1-10% of their intrinsic fluxes. The typical 
2-10 ke V reflection-dominated spectrum is flat and characterized 
by a prominent (EW i 0.6-1 keV) Fe Ka emission line at 6.4 
keV. 
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According to the X-ray background synthesis models (e.g. 
Gilli et al. 2007) and direct estimates from local samples of 
galaxies (RisaUti et al. 1999; Cappi et al. 2006; MaHzia et 
al. 2009; Burlon et al. 2011), CT AGNs should account for 
~ 20-30% of the entire AGN population. Nonetheless, only 
a few dozen of them (mostly local, low-luminosity Seyfert 
galaxies) have been unambiguously discovered and studied so 
far (Comastri 2004; Fukazawa et al. 2011; Della Ceca et al. 
2008 and references therein). Furthermore, ultra-deep X-ray 
exposures have provided clear evidence of powerful CT quasars 
at high redshifts (e.g. Feruglio et al. 2011; Comastri et al. 
2011; Gilli et al. 2011; Fiore et al. 2011), whose existence 
has been debated for many years (e.g. Halpern et al. 1999). 
Enlarging the sample of well-studied CT AGNs in the X-ray 
regime is therefore useful to characterize their properties more 
accurately, confirm distinctive spectral features, and strengthen 
the relations between optical/IR and X-ray observable quantities 
of this elusive class of X-ray sources. 



In this paper, we study the X-ray spectral properties of ESO 
138-Gl, a Seyfert 2 galaxy at z = 0.0091. This AGN is hosted 
in an E/SO galaxy that has a peculiar morphology of a bulge sur- 
rounded by a ring, within which there is a hint of a bar (AUoin 
et al. 1992; Ferruit et al. 2000). Furthermore, ESO 138-Gl 
exhibits a compact nucleus and a bright asymmetric, wedge- 
shaped circumnuclear zone of diffuse light resembling an ion- 
ization cone from the AGN (Munoz-Marin et al. 2007). All these 
features make ESO 138-Gl a unusual and particularly inter- 
esting object. On the basis of ASCA data, Collinge & Brandt 
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(2000) found a hard spectrum (F 0.5) and a strong Fe Ka 
emission line, implying that this Seyfert 2 galaxy is Compton- 
thick/reflection-dominated. However, statistically acceptable fits 
can also be achieved using a partial-covering absorption column 
of Nh = 2 X 10^^ cm"^. Here, we present and discuss the results 
of our analysis of two archival XMM-Newton spectroscopic ob- 
servations of ESO 138-Gl. 



2. XMM-Newton observations and data reduction 

ESO 138-Gl was observed by XMM-Newton (Jansen et al. 
2001) on 2007 February 16 (Obs. ID: 0405380201) and 2007 
March 25 (Obs. ID: 0405380901). Both observations were per- 
formed with the EPIC cameras operating in full-frame mode. 

Data were reduced with SAS 11.0 (Gabriel et al. 2004) us- 
ing standard procedures. We selected X-ray events correspond- 
ing to patterns 0-4(0-12) for the PN(MOS) camera. The event 
lists were filtered to ignore periods of high background flaring 
according to the method presented in Piconcelli et al. (2004) 
based on the cumulative distribution function of background 
light-curve count-rates. For the FN detector, source counts were 
extracted from a circular region with a radius of 30 arcsec (14 
arcsec for the 901 observation since it was affected by a higher 
background level). Background counts were estimated from a 
50 arcsec -radius region on the same chip for both datasets. After 
screening the final net exposure times are 15.5 and 11.4 ks for 
the 201 and 901 observation, respectively. The redistribution ma- 
trix files and ancillary response files were created using the SAS 
task RMFGEN and ARFGEN, respectively. During the XMM- 
Newton observations, ESO 138-Gl did not display any varia- 
tion in either X-ray flux and the spectral shape, thus we coadded 
the FN spectra to increase the signal-to-noise ratio and created a 
combined response matrix. Source and background spectra were 
summed using the FTOOL task mathphc^ 

The combined FN spectrum was rebinned so that each en- 
ergy bin contained at least 25 counts to allow us to use the 
minimization technique in spectral fitting. 

The reduction of MOS spectra was carried out following the 
same procedure applied to the FN datasets. We adopted an ex- 
traction radius of 31 arcsec for the source spectra, while the 
background counts were extracted from neighbouring source- 
free circular regions with a radius of 50 arcsec. 

In this paper, we focus on the analysis of FN spectroscopic 
data only because of the higher sensitivity of this detector over 
the broad 0.5-10 keV range than both MOS cameras (even when 
they are combined), and above 5 keV in particular Nevertheless, 
we checked that consistent results were obtained including the 
MOS data in our analysis. Finally, no useful high resolution RGS 
data (den Herder et al. 2001) were obtained from the present 
XMM-Newton observations of ESO 138-Gl because of their 
short exposure times. 

3. X-ray spectral analysis 

We analyzed the spectra with the XSFEC vl2.5 package 
(Arnaud 1996). We limited our analysis to the 0.5-10 keV en- 
ergy range, where the accuracy of the EPIC calibrations is max- 
imal. All models presented in this paper include absorption 
caused by the line-of-sight Galactic column density of Nh =1.3 
X 10^1 cm-2 (Kalberla et al. 2005). 



' Spectral files were combined following the method outlined in 
\http://heasarc. gsfc. nasa. gov/docs/asca/bgd_scale/bgd.scale.htmL 



We provide best-fit parameter values in the source frame, un- 
less otherwise specified. The quoted errors in the model param- 
eters correspond to a 90% confidence level for one interesting 
parameter (Ax^ = 2.71; Avni 1976). In this work, we assume a 
ACDM cosmology with Hq = 70 km s"' Mpc"' and Qa = 0.73 
(e.g., Larson et al. 201 1). 

3. 1 . Hard X-ray band and Fe Ka lines 

The 0.5-10 keV EPIC FN spectrum of ESO 138-Gl is shown 
in Fig. la. Owing to its complexity, we first studied the X-ray 
spectrum in the 3.5-10 keV band to gain insight into the origin 
of the X-ray emission in this energy range. This emission should 
be relatively unaffected by the extended photoionized/scattered 
"soft excess" component frequently observed in obscured AGNs 
(Turner et al. 1997; Sambruna et al. 2001; Bianchi et al. 2007), 
which we consider in Sect. 3.2. As expected for a highly ab- 
sorbed Seyfert 2 galaxy, a fit to the FN spectrum with a simple 
power-law in the 3.5-10 keV band (excluding the 6-7.5 keV Fe 
K energy range) reveals a very flat spectral shape with a pho- 
ton index F ~ 0.35. We tested two alternative scenarios for the 
hardX-ray spectrum of ESO 138-Gl: (i) a reflection-dominated 
(RD) one, where the X-ray primary continuum from the AGN 
is totally depressed in the EPIC band, as expected in case of a 
CT absorber ; and (ii) a transmission-dominated (TD) scenario, 
where the AGN emission is absorbed by a Compton-thin (i.e., 
Nh < 10^'* cm"-) obscuring screen that only moderately sup- 
presses the X-ray continuum in the 3.5-10 keV range. Moreover, 
we explore a scenario wherein both reflection and transmission 
components are present. 

In Fig. lb, we show the hard X-ray FN spectrum fitted by 
the best-fit RD model, which consists of a cold reflection com- 
ponent, modeled by PEXRAV in XSFEC. In this modelling, we 
assumed that F = 1.8, that the reflection scaling factor R = -I 
(defined as R - Q./2n, where Q is the solid angle subtended 
by the reflector for isotropic incident emission; Magdziarz & 
Zdziarski 1995), that the metal abundances of the reflector were 
fixed to their solar values and that the inclination angle was fixed 
to 65 deg. The RD model also included a strong fluorescent Fe 
Ka emission line close to 6.4 keV, plus three narrow Gaussian 
lines at 6.63, 7.058, and 7.54 keV associated with He-like Fe 
Ka, Fe K/3, and Ni Ka emission, respectively. The addition of 
each Fe and Ni Ka line produces an improvement in the result- 
ing fit statistic significant at >99% confidence level (94% in the 
Fe line case) according to an F-test. We fitted the Compton 
shoulder (CS) to the Fe Ka line by adding a Gaussian emis- 
sion line at 6.3 keV with o" =40 eV and a normalization fixed to 
20% of the neutral Fe Ka line as expected for CT material (e.g.. 
Matt et al. 2002; Yaqoob & Murphy 201 1). The CS redwards of 
6.4 keV is emitted from neutral material as a result of Compton 
down-scattering of the Fe line photons in the same medium. The 
addition of the CS component reduces the fit statistic of Ax^ 
- 6. The RD model gives a very good fit to the data with ax^iv) 
= 78(75) (see Table 1). 

The TD model yields a reasonable fit to the XMM- 
Newton data with a final x^(^) - 83(73). The spectrum is de- 
scribed by an absorbed continuum power-law with F = 1.9+0.4 
and Nh = 2.4+^-^ x 10^^ cm"^ (see Table 1). As for the model 
RD, four narrow Gaussian lines (i.e., neutral Fe Ka, He-like Fe 
Ka, Fe KJ3, and Ni Ka) are included in the fit. When we add a 
CS component to the fit that a normalization equal to 10% of the 
Fe Ka line core (as expected for a column density of ~ 2.5 x 
10^-^ cm"^; Yaqoob & Murphy 201 1), we achieve a Ax^ = 4. 
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Energy (keV) Energy (keV) 

Fig. 1. (a)-Left: XMM-Newton PN spectrum resulting from the merging of the 2007 February and March observations of ESO 
138-Gl. (b)-Right: The RD model fitted to the 3.5-10 keV spectrum of ESO 138-Gl. This model consists of a Compton reflection 
continuum component and five Gaussian lines accounting for the neutral Fe Ko' line, its Compton shoulder, and the He-like Fe Ka, 
Fe K/3, and Ni Ka emission lines, respectively. The bottom panel shows the data-to-model ratios. 



We investigated whether a clumpy absorber might be present 
along our line of sight by adding an absorbed power-law compo- 
nent to the TD model. The two power laws had the same photon 
index, thus mimicking a scenario wherein the obscuring screen 
along our line of sight is not uniform (e.g., Elitzur 2008; Ramos- 
Almeida et al. 2009). The fit resulted inx^iv) = 75(71), with an 
unusual steep continuum (T - 2.8+0.4) seen through two distinct 
layers of absorption of cold gas with column densities of ~ 7.7 
and ~ 1.5 X 10^^* cm"^, respectively. However, the spectral pa- 
rameters of the dual absorber model were largely unconstrained. 
We then fit the data with the continuum slope fixed to F = 1.9 
(as found for the TD model). In this case, the column density 
of one absorption component was found to be Nh ~ 3.5 x 10^-* 
cm"^, while the Nh value for the other component was found to 
be consistent with zero. The ratio of the normalization of two 
power-law components is ~9. 

The hard continuum can also be reproduced with a more 
complex model where reflection and transmission are both 
present, which provides a fit Wilhx^iy) - 75(72) when both Nh 
and F are allowed to vary. The ratio of the normalization of the 
reflection to that of the transmitted component is ^ 6, implying 
that the reflection continuum can easily account for most of the 
hard X-ray emission in ESO 138-Gl. However, we note that 
the quality of our data, along with the lack of spectral coverage 
above 10 keV, does not allow to properly constrain the relative 
contribution of both components. Fixing the value of the column 
density of the absorber to Nh = 1.6 x lO^'* cm"^ (= cr,"') yields a 
flattening of the photon index to F ~ 1.85, but does not improve 
the quality of the fit statistics. 

As performed by Guainazzi et al. (2005) and LaMassa et al. 
(201 1) for their samples of heavily obscured Seyfert 2 galaxies, 
we performed a "local" fit to the spectrum between 5.5 and 7 
keV using Cash statistics (Cash 1979), which allows the use of 
unbinned data without any loss of spectral information, to ob- 
tain a more robust determination of the neutral Fe Ka line pa- 
rameters. The continuum underlying the emission line was mod- 
eled with a power law, with the photon index fixed to the best-fit 
value of F = 0.35. The line energy centroid is at 6.432+0.015 
keV, with an upper limit to the line width of cr < 60 eV. The EW 
measured with respect to the underlying continuum is 790±70 
eV (see Table 2), including the contribution from the CS (which 
originates from Compton backscattering of line photons). By fit- 
ting the MOS data, the properties derived for the Fe Ka line are 
fully consistent with those measured by the PN spectrum. 



Table 1. Best-fit spectral parameters for the hard X-ray (3.5-10 
keV) continuum of ESO 138-Gl. See Sect.Ofor details. 



Parameter 


Model RDf 


Model TD f 


P 


1.8* 


1 Q+0.4 


Norm* 
N^ 

fd 

' 2-10 

;r'(dof)^ 


o '-)+0.4 

2.3 
4.3 
1.04(75) 


16^^^^ 

2.4i! 

2.1 
9.6 
1.14(73) 



* Spectral model: RD (re/Zection-dominated), TD (transmission- 
dominated); " photon index of primary continuum; * normalization of 
the reflection (RD) or absorbed power-law component (TD) (10^^^ pho- 
tons keV"' cm"2 s"' at 1 keV); ^ column density of the absorber (10^^ 
cm-2); 2-10 keV flux (lO^'^ erg cm'^ s"'); " 2-10 keV luminosity 
(lO"*' erg s"'); ^ reduced and number of degrees of freedom. ' Fixed 
value (see text for details). 

This value of the EW of the 6.4 keV Fe Ko- emission line un- 
ambiguously rules out the TD scenario (as well as the Compton- 
thin dual absorber model) because it is too large for a contin- 
uum obscured by the Compton-thin column density of a few 10^^ 
cm"^ derived from the TD fit to the data (e.g., Ghisellini, Haardt 
& Matt 2004; Guainazzi et al. 2005; Fukazawa et al. 201 1). The 
prominent Fe Ko- line observed in ESO 138-Gl lends strong 
support to the view that the X-ray primary continuum from the 
AGN is totally depressed in the EPIC band, as expected for the 
presence of a CT absorber, and the 2-10 keV emission is due to 
the reprocessing of the AGN light by material surrounding the 
nucleus. Consequently, we adopt hereafter the RD model as the 
best-fit description of the hard X-ray spectrum of ESO 138-Gl. 

3.2. Broad-band spectroscopy 

As shown in Fig. la, the soft X-ray portion of the XMM-Newton 
spectrum of ESO 138-Gl displays several emission features, 
in addition to a smooth excess below ~ 3 keV deviating from 
the flat shape observed for the hard X-ray band. As commonly 
done for low-resolution spectra of heavily obscured AGNs (e.g.. 
Pounds & Vaughan 2006; Piconcelli et al. 2007), we fitted this 
"soft excess" with a phenomenological model consisting of the 
RD model, an unabsorbed power-law component, and a se- 
quence of narrow Gaussian lines to account for the emission 
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Table 2. Best-fit spectral parameters for the emission lines de- 
tected in the hard X-ray band. See Sect. B.ll for details. 



Line Id. Energy EW Intensity 

(1) (2) (3) (4) 

FeKa* 6.432+^-^;;^ 790+70 2.34±0.25 x 10"' 

Fe XXV 6.631"' ^ i25+ll 4.5+J° x IQ-'' 

Fe K/3 7.058'" 80+*» 2.3!|-^ x lO"* 

NiKo- 7.54+ 150+90 2.9+1.8 x 10-* 



The columns give the following information: (1) line identification; (2) 
energy (ke'V); (3) E'W (e'V); and (4) intensity (photons cm^^ s"') of the 
line. ' Properties of the Fe Ka emission line are derived by a 'local' fit 
of the unbinned spectrum in the 5.5-7 keV energy range, using the Cash 
statistics. The EW value includes the contribution of the CS. ' Fixed 
value. 

features observed in the ~ 0.5-2.5 keV range. The power law 
was found to be steep (Fsoft = 3.0+^^) and appears to describe 
a mixture of emissions caused by the electron-scattered fraction 
of the primary continuum, thermal plasma, and blends of un- 
resolved emission lines and radiative recombination continua, 
as observed in most X-ray obscured AGNs (Turner et al. 1997; 
Kinkhabwala et al. 2002). Table 3 lists the best-fit model param- 
eters for the emission lines detected in the soft X-ray band. These 
emission lines are associated with atomic transitions commonly 
detected in the soft X-ray spectra of obscured AGNs. However, 
owing to the limitations imposed by spectral resolution of EPIC 
and potential line blending with adjacent emission lines, both 
the line identification and parameters should be interpreted with 
care. 

This model provides a good fit to the X-ray spectrum of ESO 
138-Gl in the 0.5-10 keV range (Fig. 2), with an associated 
xHv) = 156(133). 

An alternative explanation of the "soft excess" in some ob- 
scured AGNs is emission related to star-forming activity. Here 
we used a combination of two thermal-emission components 
(MEKAL model in XSPEC) to fit the soft X-ray portion of the 
spectrum of ESO 138-Gl. We estimated a metallicity of Z/Z© - 
0.3±0.1 and a temperature of kTi « 0.7 and kT2 < 0.1 keV, re- 
spectively. However, this model yields a poor fit with^J(dof) - 
1.6(152), leaving significant positive residuals in the 1-2.5 keV 
region and, thus, disfavoring a scenario wherein the soft excess 
in ESO 138-Gl originates in thermal emission from a starburst. 
This is an expected result because a substantial thermal contribu- 
tion to the soft X-ray emission is typically found only in Seyfert 
2/starburst composite galaxies (e.g., Levenson et al. 2001a,b), 
while the host galaxy of ESO 138-Gl does not exhibit evidence 
of significant star-forming activity at optical wavelengths (AUoin 
et al. 1992; Ferruit et al. 2000). 

4. Fluxes and luminosities 

Once our broad-band best-fit model is assumed, we measure a 
hard X-ray flux F2-10 = 2.3 x 10"'^ erg cm"^ s"' and a soft X- 
ray flux ^0.5-2 = 3.1 X lO"'-' erg cm"^ s"'. After correction for 
Galactic absorption, they correspond to an observed luminosity 
L2-10 = 4.3 X lO'*' erg s"' and Lo.5-2 = 8.8 x 10"*" erg s"' in the 
hard and soft band, respectively. 

The T ratio of the observed 2-10 keV to the extinction- 
corrected [OIII] fluxes can be used as a proxy for the amount 
of obscuration of the X-ray primary continuum, since the [OIII] 



Table 3. Best-fit spectral parameters for the emission lines de- 
tected in the soft X-ray band. See Sect. 3.2 for details. 



Energy 


Intensity 


Identification 


(1) 


(2) 


(3) 


0.57+0.02 


5.8^3 Q X lO""" 


OVII LyQ- 


0.72+0.01 


3.0+1.3 X 10 


OVII RRC, Fe XVII 3s-2p 


U.o4+U.Uz 


1 ^_i_n T ^ 1 n— 5 
Z.j+y). / X lU 


re AVli-AVili ja-zp, UVill KKU 


0.93+0.02 


3.3+0.9 X 10-5 


NelX Hea 


1.02"f 


1.2+1;^ X 10~5 


NeX Lya, Fe XXI 3d-2p 


1.14+0.03 


9.3+5.8 X 10-« 


Fe XXIIl-XXIV L 


1 91+0.04 
^•^^-0.02 


1.1+0.5 X IQ-^ 


NeX Ly/3 


1.34+0.02 


1.2+g 2 X 10-5 


MgXI Hea 


1.50+0.02 


6.8+^2 X 10-" 


MgXII Lya 


1.82+0.02 


7.6+1'^ X 10-" 


SXIII Hea 


7 04+0.08 


2.0+|-^ X I0-" 


Si XXIV-XVI Lyo- 


7 4^+0.02 


6.2+1.8 X 10-" 


SXV Hea 



The columns give the following information: (1) energy of the line 
(keV); (2) intensity of the line (photons cm-^ s-'); (3) likely identifica- 
tion (e.g., Brinkman et al. 2002; Guainazzi & Bianchi 2007). + indicates 
that the energy has been fixed to the laboratory energy of NeX Lya due 
to the likely superposition of many unresolved Fe lines. 



line is an isotropic indicator of AGN power as it is produced 
in the NLR. In particular. X-ray sources for which T < 1 are 
typically associated with a CT absorber and a 2-10 keV RD 
spectrum (Bassani et al. 1999; Akylas & Georgantopoulos 2009; 
Gonzalez-Martin et al. 2009). The extinction-corrected [OIII] 
flux of ESO 138-Gl is 2.7 x 10-'^ erg cm-^ s-', being de- 
rived from the observed value of 9.72 x 10-'^ erg cm-^ s-' 
reported in Schmitt & Storchi-Bergmann (1995) after the cor- 
rection (e.g., Lamastra et al. 2009 and references therein) for a 
Balmer decrement Ha/H/i = 4.25 (Alloin et al. 1992). This gives 
T = 0.85, providing further support to the existence of a CT 
screen that suppresses the observed 2-10 keV flux. Furthermore, 
in the framework of a TD scenario for the hard X-ray spectrum 
of ESO 138-Gl, we note that the measured value of T is in- 
consistent with a Compton-thin column density of ~ 2.5 x 10^^ 
cm-^. In this case, the absorption-corrected X-ray flux would 
be F2-10 = 5.1 x 10-'^ erg cm-- s-' and, therefore, T is ~ 
1.9. However, as found by many studies (Maiolino et al. 1998; 
Panessa et al. 2006; Lamastra et al. 2009), unabsorbed Seyfert 
1 and absorption-corrected Compton-thin Seyfert 2 galaxies in- 
deed typically show T > 10. 

5. Discussion 

We have presented the analysis of the combined spectroscopic 
data from two XMM-Newton observations of the Seyfert 2 
galaxy ESO 138-Gl performed in 2007. The 0.5-10 keV spec- 
trum of this source is typical of a heavily obscured AGN: a 
very flat hard X-ray continuum with a strong Fe Ka emission 
line at 6.4 keV and a "soft excess" component characterized 
by many emission lines from highly ionized metals. In partic- 
ular, two main pieces of evidence lead us to consider the RD 
scenario as the most physically plausible description for the 
XMM-Newton data of ESO 138-Gl, thus confirm the claim by 
Coflinge & Brandt (2000) based on 1997 ASCA data for the ex- 
istence of a CT absorber covering the primary X-ray continuum 
in this Seyfert galaxy. The first piece of evidence is the value of 
EW ~ 800 eV for the Fe Kcc line around 6.4 keV. Such a large 
value cannot be explained if the underlying continuum is trans- 
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Energy (keV) Energy (keV) 

Fig. 2. (a)-Le/f; The 0.5-lOkeV spectrum of ESO 138-Gl when the best-fit model is applied. The lower panel shows the deviations 
of the observed data from the model in units of standard deviation. (h)-Right: Best-fit model for the broad band X-ray spectrum. This 
model (solid line) consists of a pure reflection component (dash-dotted line) plus an additional unobscured power-law component 
accounting for the soft X-ray scattered/leaked emission (dashed line). It also includes the hard and soft X-ray Gaussian emission 
lines listed in Tables 2 and 3, respectively. 



mitted through a Compton-thin absorber with Nh ~ 2.5 x 10 
cm"^, as inferred by the TD fitting model. A strong Fe Ka emis- 
sion line with an EW i 0.7-1 keV is indeed interpreted as the 
unambiguous signature of a hard X-ray emission produced by re- 
processing the intrinsic X-ray continuum and then scattering this 
radiation along the line of sight (Ghisellini, Haardt & Matt 2004; 
Guainazzi et al. 2005). In this RD scenario, an obscuring screen 
with Nh i 10^"^ cm"^ is believed to completely block the X-ray 
primary source from direct view. It has been generally assumed 
that the absorber and the reflector consist of the same material, 
i.e. the torus, homogeneously fills a doughnut-shaped region lo- 
cated between the BLR and NLR (Matt et al. 1996; Molendi et 
al. 2003). However, the discovery of fast transitions from a CT to 
a Compton-thin spectral state (or vice-versa) have been observed 
in a handful of Seyfert galaxies (Risaliti et al. 2005; Bianchi et al. 
2009) and interpreted in terms of absorbing clouds crossing our 
line of sight and located very close to the X-ray source (Elitzur 
2008). This implies that a line-of-sight clumpy sub-pc-scale ab- 
sorber exists and, most importantly, the coexistence in the same 
source of multiple absorbing and/or reflecting components dis- 
tributed across a range of distances from fractions of pc up to 
several pc from the SMBH. However, there has so far been no 
evidence in ESO 138-Gl of large variations in the CT absorber 
Future X-ray observations of this source might detect RD to TD 
spectral state transitions, hence reveal both the presence of an in- 
homogeneous obscuring gas and the intrinsic luminosity of the 
X-ray continuum. Finally, the larger EW value of 1.7+0.4 keV 
for the line at 6.4 keV reported by Collinge & Brandt (2000) can 
be explained in terms of a blend of the lines associated with Fe 
Ka emission, caused by the poorer spectral resolution of ASCA 
detector than that provided by the PN camera. 

Secondly, the ratio T = 0.85(1.9) of the observed(intrinsic) 
2-10 keV flux to the reddening-corrected [Olll] emission line 
flux is consistent with a RD X-ray spectrum (e.g.. Sect. 4), since 
this value does not agree with those reported for X-ray sources 
attenuated by a column density Nh < 10^"^ cm"^ (Maiolino et al. 
1988; Bassani et al. 1999; Heckman et al. 2005; De Rosa et al. 
2008). Furthermore, the [OIII] luminosity can be used to infer 
the intrinsic 2-10 keV luminosity of the AGN in ESO 138-Gl. 
Using the Eq. (3) in Lamastra et al. (2009), we derive a L2-10 ~ 
6 X 10"^^ erg s"', i.e. well within the typical Seyfert range (and 
a factor of six larger than the value inferred by assuming the 
Compton-thin/TD scenario). The ratio of the value of the intrin- 



sic 2-10 keV luminosity estimated by the [OIII] to the observed 
value measured from the RD model is =s 14, which is typical of 
CT AGN (Comastri 2004; Turner et al. 1997). A value of L2-10 
X 10*^ erg s"' for ESO 138-Gl can also be inferred from its 
mid-IR luminosity (i.e., Lmir = 3.46 x lO"^-' erg s"'), as shown 
by Gandhi et al. (2009). We note that the X-ray luminosity cal- 
culated with the TD model, i.e. 9.6 x lO'*' erg s ' (see Table 
1), is about an order of magnitude lower than this MIR-based 
estimate. 

The estimate of the intrinsic hard X-ray luminosity of ESO 
138-Gl allows us to derive some information about the SMBH 
mass of this object. For a [OIII]-based La-io ~ 6 x 10"*^ erg s"\ 
Vasudevan & Fabian (2007) indicate a bolometric correction of 
10. Accordingly, we infer an Eddington luminosity LEdd = 6 x 
lO**"* erg s"' by assuming a ratio L/LEdd = 0.1, which is typical 
of Seyfert galaxies in the local Universe (Woo & Urry 2002). 
This value of the Eddington luminosity implies a black hole 
mass of 4.6 x 10* Mq. It is well-known that near-infrared bulge 
luminosities and SMBH masses show a tight correlation (e.g., 
Marconi & Hunt 2003). The /T-band magnitude of the bulge 
of ESO 138-Gl is K^IO.71 mag (Peng et al. 2006), which 
translates into a /iT-band luminosity of Lx^Buige - 6.6 x 10"^ 
Lq. Accordingly, we find that ESO 138-Gl does not lie on the 
MBH-LK,Buige Correlation: the estimated mass of the SMBH is 
about a factor of ten lower than expected. This implies that ESO 
138-Gl is a peculiar early-type galaxy with an under-massive 
black hole, or an intrinsic X-ray luminosity ten times higher than 
that estimated based on the [OIII] luminosity, thus providing 
further support for a CT nature of this AGN. 



Finally, useful insights into the nature of the obscuring 
screen along the line of sight to the nucleus in ESO 138-Gl 
can also be derived from the 15-150 keV flux level constraint. 
We note that ESO 138-Gl is indeed included in the catalogs 
of hard X-ray sources detected by INTEGRAL IBIS (Bird et 
al. 2010) and Swift BAT (Cusumano et al. 2010) above 15 
keV. Specifically, it is associated with the BAT source 2PBC 
J1651.9-5914, for which F15-150 = 2.3±0.5 x 10"" erg cm^^ 
s"', that is centered at a distance of ~ 3 arcmin from ESO 
138-Gl. On the basis of the RD model (i.e. assuming a broad- 
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band, pure reflection spectrum), it is possible to estimat^ a 15- 
ISOkeVfluxof 3.1!JJ-^ x 10"" erg cm^^ s"'. 

However, since the angular resolution of BAT is ~ 20 ar- 
cmin FWHM and the X-ray loud starburst/Seyfert 2 composite 
galaxy NGC 6221 (e.g. Levenson et al. 2001a) is located at 11 
arcmin from our target, the BAT flux can only be used as an up- 
per limit to the X-ray emission from ESO 138-Gl in the 15-150 
keV band (V. La Parola, private communications). The same ar- 
gument is also valid for the INTEGRAL IBIS flux reported in 
Bird et al. (2010). We were also able to extract the PN spectrum 
of NGC 6221, as it falls within the field of view of the detector 
(although close to the edge) in our observations, to help us de- 
termine the possible contribution of this source to the observed 
BAT flux. The X-ray spectrum of NGC 6221 is weU-fitted by a 
combination of a soft thermal plasma component (kl ~ 0.6 keV) 
and a Compton-thin (Nh ~ 10^' cm"^) absorbed power law with 
r = 1 .6+0. 1 , as found by Levenson et al. (2001 a) based on ASCA 
datc0. By extrapolating this model to energies > 10 keV, we esti- 
mate a value of Fis-iso = 1.1 x 10"" erg cm"^ s"' for the flux in 
the 15-150 keV band. Therefore, a sizable fraction of the BAT 
flux ascribed to ESO 138-Gl in the Cusumano et al. catalog 
might be due to the emission from NGC 6221: in particular, it is 
likely that the total 15-150 keV flux of Fi5_i5o ~ 2.3 x 10"" erg 
cm"^ s"' seen by BAT and IBIS is produced roughly equally by 
the two sources. In turn, this may indicate that the column den- 
sity of the nuclear absorber in ESO 138-Gl should be Nh i 
10^^ cm"^, otherwise any additional emission component in the 
15-150 keV band caused by the X-ray primary continuum trans- 
mitted through an absorber with Nh ~ a few 10^"^ cm"^, would 
produce a further increase in the final flux, significantly exceed- 
ing the IBIS/BAT flux values. 

In this respect, future observations of ESO 138-Gl above 
10 keV will be crucial to shed light on the properties of the 
X-ray continuum emission and the obscuring gas in the nuclear 
environment of this Seyfert 2 galaxy. We note, however, that 
given the proximity of NGC 6221, low -resolution (i.e. with 
FWHM » 10 arcmin) observations carried out with current 
X-ray observatories such as Swift BAT, INTEGRAL IBIS and 
Suzaku PIN cannot provide an unambiguous and definitive test 
of whether ESO 138-Gl contains a heavily CT (i.e. Nh » lO^'* 
cm"^) absorber. This will be achieved by next generation X-ray 
missions such as NuSTAR (Harrison et al. 2010) and ASTRO-H 
(Takahashi et al. 2010), which should to fly in 2012 and 2014, 
respectively, as they combine ~ 50-90 arcsec imaging with high 
spectral resolution in the ~ 10-60 keV energy range. 

The soft portion of the X-ray spectrum of ESO 138-Gl 
shows the typical properties of other obscured AGNs, with the 
presence of a highly structured "soft excess". High-resolution 
imaging and spectroscopy of Seyfert 2 galaxies performed with 
XMM-Newton and Chandra have allowed us to reveal that the 
soft X-ray spectrum originates from extended emission pho- 
toionized by the AGN (Kinkhabwala et al. 2002; Bianchi et al. 

^ We note, however, that any accurate estimate of the reflected emis- 
sion depends on the geometry of the reflection medium (e.g., Ikeda et 
al. 2009), which is basically unknown. Meaningful constraints on the 
geometry of the reprocessor will be possible only using the high signal- 
to-noise ratio, broad-band spectral data expected from observations per- 
formed by ASTRO-H and NuStar in the near future. 

3 The 2-10 keV flux of NGC 6221 measured by the XMM-Newton 
observation is F2-10 ~ 3 x 10"'^ erg cm"- s"', i.e. a factor of 4.5 lower 
than the value derived by previous ASCA data, thus confirming the pres- 
ence of large flux variability in this source as reported in Levenson et 
al. (2001b). 



2006; Guainazzi & Bianchi 2007). The contribution of colli- 
sionally ionized plasma associated with starburst regions does 
not appear to be significant, except in objects known to host in- 
tense star-formation activity, such as highly disturbed systems 
and very luminous IR galaxies (Netzer et al. 2005; Guainazzi et 
al. 2009; PiconcelH et al. 2010). NaiTow-band [OIII] images of 
ESO 138-Gl show a jet-like feature extending up to a; 2 kpc 
westwards from the center, produced by the ionizing nuclear ra- 
diation (Schmitt & Storchi-Bergmann 1995). Furthermore, high- 
resolution near-UV and [OIII] images of ESO 138-Gl taken 
with the ACS and WFPC2 detectors onboard HST (Munoz- 
Marin et al. 2007; Ferruit et al. 2000) have revealed a bright, 
ionization cone-shaped circumnuclear zone (~ 200 pc) of dif- 
fuse light. Close morphological correspondences between the 
soft X-ray and the [OIII] emission on a kpc scale are frequently 
observed in heavily obscured Seyfert 2 galaxy and clearly in- 
dicate a common physical origin (i.e. photoionization) for the 
two emissions (Brinkman et al. 2002; Bianchi et al. 2006 and 
references therein). Accordingly, we suggest that this may turn 
out to also be the case for the origin of the line-rich "soft ex- 
cess" observed in the XMM-Newton spectrum of ESO 138-Gl. 
However, high-resolution X-ray spectroscopy and imaging with 
deep exposures of ESO 138-Gl are needed to unambiguously 
confirm the presence of photoionized gas in this Seyfert 2 galaxy 
and resolve its morphological details. 
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